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PREFACE 


In  Hay  1979,  storm-generated  waves  caused  significant  damage  to  a  newly 
completed  diked  dredged  material  disposal  facility  located  in  Saginaw  Bay, 
Michigan.  As  a  consequence,  the  U.  S.  Army  Engineer  Waterways  Experiment  Sta¬ 
tion  (WES)  was  requested  by  the  U.  S.  Army  Engineer  District,  Detroit  (NCE), 
to  develop  a  wave-measurement  system  and,  concurrently,  a  wave  hindcast  model 
in  order  to  determine  the  frequency  if  occurrence  of  severe  wave  conditions  in 
Saginaw  Bay.  Funds  were  authorized  by  NCE  on  2  July  1979.  The  study  was  con¬ 
ducted  during  the  period  from  July  1979  to  February  1983  in  the  Wave  Processes 
Branch,  Wave  Dynamics  Division,  Hydraulics  Laboratory,  WES,  under  the  direc¬ 
tion  of  Hr.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory,  Dr.  R.  W. 
Whalln,  former  Chief  of  the  Wave  Dynamics  Division,  and  Hr.  C.  E.  Chatham, 
Acting  Chief  of  the  Wave  Dynamics  Division. 

The  study  was  conducted  under  the  direction  of  Mr.  A.  W.  Garcia,  who 
prepared  the  wave  gaging  portion  of  the  report.  Dr.  R.  E.  Jensen  developed 
the  shallow- water  wave  hindcast  model  and  prepared  the  theoretical,  compari¬ 
son,  and  design  wave  condition  sections  of  the  report.  Also,  M.  B.  Habeeb 
and  B.  F.  Vavra  coordinated  all  activities  associated  with  the  final  publica¬ 
tion  of  this  report. 
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Commanders  and  Directors  of  WES  during  the  conduct  of  the  study  and  the 
|  preparation  and  publication  of  this  report  were  COL  John  L.  Cannon,  CE, 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  can  be  converted  to  metric  (SI)  units 


as  follows: 

Multiply 

By 

To  Obtain 

acres 

4046. 856 

square  metres 

cubic  yards 

0.7645549 

cubic  metres 

feet 

0.3048 

metres 

feet  per  second 
per  second 

0.3048 

metres  per  second 
per  second 

knots  (international) 

0.514444 

metres  per  second 

miles  per  hour 
(U.  S.  statute) 

1.609344 

kilometres  per  hour 

miles  (U.  S.  nautical) 

1.852 

kilometres 

miles  (U.  S.  statute) 

1.609344 

kilometres 

LAKE  HURON 


WAVE  DATA  ACQUISITION  AND  HINDCAST  FOR 
SAGINAW  BAY.  MICHIGAN 


PART  I:  INTRODUCTION 

1.  The  Saginaw  River  is  a  navigable  waterway  that  flows  northward  into 
the  head  of  Saginaw  Bay,  Michigan.  The  areas  along  the  river  and  its  tribu¬ 
taries,  the  Flint,  Cass,  Tittabawassee ,  and  Shiawassee  Rivers,  are  heavily 
industrialized;  and  for  years,  industrial  and  sewage  treatment  wastes  have 
been  discharged  into  these  rivers  resulting  in  contamination  of  the  bottom 
sediments.  At  the  request  of  the  Governor  of  Michigan,  maintenance  dredging 
of  the  lower  river  and  Saginaw  Bay  navigation  channel  was  discontinued  in 
1970.  Dredging  was  suspended  pending  the  selection  of  a  site  or  construction 
of  a  containment  facility  suitable  for  disposal  of  the  contaminated  dredged 
material.  In  conjunction  with  the  Bay  County  Michigan  Board  of  Commissioners, 
the  U.  S.  Army  Engineer  District,  Detroit  (NCE) ,  selected  a  site  to  the  east 
of  the  mouth  of  the  Saginaw  River  as  the  location  for  construction  of  a  diked 
dredged  material  disposal  facility.  Construction  was  begun  in  1977  and  com¬ 
pleted  in  the  fall  of  1978  at  a  cost  of  about  14  million  dollars.  The  facil¬ 
ity  covers  an  area  of  about  285  acres*  in  water  depths  of  10  to  14  ft  and  is 
protected  by  a  rubble-mound  breakwater.  The  design  capacity  is  approximately 
10  million  cubic  yards  with  an  initial  fill  of  4  million  cubic  yards  and 
600,000  cubic  yards  yearly  for  10  years,  the  anticipated  active  life  of  the 
structure.  Upon  filling  of  the  facility,  the  resulting  island  will  revert  to 
Bay  County  Authorities,  probably  for  recreational  use.  Figure  1  shows  the 
location  and  scale  of  the  structure.  The  cross  dike  was  necessary  for  two 
reasons.  First,  construction  of  the  structure  required  two  building  seasons. 

A  partially  completed  structure  would  have  been  vulnerable  to  wave  and  ice 
damage  during  the  fall  and  winter  months.  The  cross  dike  served  to  close  the 
partially  completed  structure  tftiile  construction  was  suspended  during  the 
winter.  Second,  the  dike  served  to  support  the  dredge  discharge  pipe  in  order 
to  assure  uniform  filling  of  the  facility. 


* 


! 


A  table  of  factors  for  converting  U.  S.  customary  units  of  measurements  to 
metric  (SI)  units  is  presented  on  page  3. 


PART  II:  STORM  OF  MAY  1979 


2.  During  the  period  23  to  25  May  1979,  a  storm  accompanied  by  sus¬ 
tained  north  to  northeast  winds  of  30  to  35  mph  produced  waves  that  caused 
significant  damage  to  the  facility.  The  damage  was  confined  to  the  northeast 
sector  of  the  structure  but  the  cost  of  repair  and  rehabilitation  was  in  ex¬ 
cess  of  one  million  dollars.  The  damage  was  primarily  caused  by  waves  over¬ 
topping  the  structure.  Figure  2  is  an  aerial  photograph  taken  on  25  May-  1979 
showing  the  structure  being  overtopped  during  the  storm.  Note  the  breach  in 
the  crown  of  the  structure  that  had  occurred  prior  to  the  photograph  being 
taken.  Figure  3  is  a  photograph  taken  on  29  May  1979,  4  days  after  the  storm. 
Roughly,  the  same  section  of  the  structure  is  shown;  note  the  damage  to  the 
crown  of  the  structure. 

3.  There  was  conjecture  following  damage  to  the  structure  that  although 
designed  in  accordance  with  guidelines  of  the  Shore  Protection  Manual  (U.  S. 
Army  CERC  1977)  for  the  known  weather  patterns  of  the  area,  the  wave  climate 
in  Saginaw  Bay  may  have  been  more  severe  than  anticipated. 

4.  As  a  consequence  of  damage  to  the  structure,  the  NCE  requested  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  to  design  and  deploy  a 
wave-measurement  system  in  Saginaw  Bay  and  concurrently  to  develop  the  capa¬ 
bility  of  hindcasting  the  wave  climate  within  the  bay.  The  wave-measurement 
system  was  initially  deployed  in  the  fall  of  1979.  It  was  configured  to  pro¬ 
vide  data  at  two  locations  within  the  bay  and  one  location  in  Lake  Huron  just 
outside  the  bay  entrance  (Figure  1) .  The  locations  referenced  to  Figure  1 
were  chosen  to  provide  data  as  follows: 

a_.  Gage  Site  1.  The  wave  climate  in  the  vicinity  of  the  disposal 
facility  primarily  during  north  and  northeast  winds. 

_b.  Gage  Site  2,  The  wave  climate  bayward  of  the  shoal  near  the 
entrance  to  the  bay  during  north  and  northeast  winds.  Also,  by 
comparing  the  wave  data  at  location  2  with  data  acquired  at  lo¬ 
cation  1,  the  wave  generation  (growth  versus  fetch  and  time) 
within  the  bay  could  be  determined . 

£.  Gage  Site  3.  The  wave  climate  just  to  the  outside  of  the  shoal 
near  the  entrance  to  the  bay  during  north  and  northeast  winds. 

By  comparing  the  wave  data  at  location  3  with  data  acquired  at 
location  2,  the  amount  of  wave  energy  being  transmitted  across 
the  shoal  from  Lake  Huron  to  Saginaw  Bay  can  be  determined. 

5.  During  the  fall  1979,  spring  1980,  and  fall  1981  deployments,  the 
wave  gages  were  set  to  acquire  17  min  of  wave  data  every  3  hr  at  a  sampling 


rate  of  2  Hz.  Examination  of  data  obtained  during  minor  storms  occurring  in 
these  deployment  periods  revealed  there  could  be  substantial  changes  in  the 
wave  climate  between  3-hr  sampling  periods.  For  the  spring  1981  deployment, 
it  was  decided  to  change  the  measurement  format  to  acquire  17  min  of  wave  data 
every  2  hr  at  a  sampling  rate  of  1  Hz.  The  modified  format  still  provided 
adequate  definition  of  3-sec  and  longer  period  waves  while  permitting  a  reduc¬ 
tion  in  the  interval  between  records  from  3  to  2  hr,  thus  providing  a  more 
definitive  record  of  the  temporal  changes  in  the  wave  climate. 

6.  It  was  assumed  at  the  inception  of  the  study  that  little  of  the  wave 
energy  propagating  toward  Saginaw  Bay  from  Lake  Huron  would  be  transmitted 
across  the  shoal.  There  was,  however,  no  observational  basis  for  this  assump¬ 
tion.  If  this  assumption  proved  correct,  Saginaw  Bay  could  be  treated  as  a 
closed  basin  for  the  purpose  of  hindcast  modeling.  This  would  eliminate  the 
need  to  hindcast  wave  conditions  in  Lake  Huron  concurrent  with  those  in  Sagi¬ 
naw  Bay  and  greatly  simplify  the  problems  involved  in  the  modeling  effort. 
Figures  4  and  5  show  spectra  obtained  at  locations  2  and  3,  respectively,  dur¬ 
ing  the  May  1981  storm.  Examination  of  these  spectra  shows  that  very  little 
of  the  wave  energy  in  Lake  Huron  is  propagating  into  Saginaw  Bay.  Indeed,  the 
position  of  the  peak  of  the  spectrum  obtained  at  location  2  implies  the  waves 
there  to  be  almost  totally  locally  generated  under  conditions  of  north  and 
northeast  winds,  the  most  severe  case  for  the  assumption  to  be  tested.  Fur¬ 
ther  examination  of  the  time-history  of  significant  height  (measure  of  the 
total  energy  at  the  gage)  results  shown  in  Figure  6  reveals  that  initial 
growth  rates  at  all  three  gage  locations  are  uniform.  Wave  conditions  at  gage 
sites  3  and  2  become  saturated  while  differences  in  the  significant  height  re¬ 
sults  vary  as  much  as  3.5  ft.  Data  in  Figures  4,  5,  and  6  clearly  show  that 
wave  conditions  found  in  Saginaw  Bay  will  be  caused  by  wave  growth  within  the 
bay. 

7.  Comparison  of  spectra  obtained  at  locations  1  and  2  allows  calcula¬ 
tion  of  the  wave  growth  as  functions  of  fetch  and  time.  The  fetch  between 
these  locations  is  about  16  statute  miles. 
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Figure  5.  Wave  spectrum  obtained  at  gage  site  3  during  May  1981  storm 


PART  III:  WAVE-MEASUREMENT  SYS  TO! 


8.  The  wave  data  acquisition  system  consisted  of  three  completely  sub¬ 
merged  quadripod  platforms,  one  at  each  of  the  locations  shown  in  Figure  1. 

On  each  platform  was  mounted  a  commercially  obtainable  SEA  DATA  Model  635  wave 
gage.  These  are  self-contained,  pressure  sensing  instruments  which  record  on 
standard  Phillips  type  magnetic  tape  cassettes.  This  type  of  instrument  was 
selected  because  the  preferred  locations  are  too  far  from  shore  to  economi¬ 
cally  permit  cabled  sensors,  and  storms  that  tend  to  cause  damage  to  the  dis¬ 
posal  facility  are  commonly  accompanied  by  lightning  which  can  interfere  with 
or  disable  a  telemetering  type  Instrument.  Moreover,  Saginaw  Bay  is  a  popular 
boating  and  fishing  area  with  the  attendant  boat  traffic.  Because  of  this 
traffic,  it  was  felt  that  surface-following  wave  buoys  would  be  vulnerable  to 
vandalism  or  theft. 

9.  The  quadripod  platforms  were  positioned  using  electronic  triangula- 

tlon.  It  was  originally  thought  that  acoustic  releases  and  beacons  would  be 
necessary  for  locating  the  platforms  in  order  to  service  the  gages.  However, 

the  electronic  positioning  proved  to  be  repeatable  to  within  about  10  m  of  the 

positions  of  the  platforms,  thereby  eliminating  the  need  for  the  releases  and 
beacons.  Small  witness  buoys  were  used  to  mark  the  location  of  the  platforms 
but  these  were  expendable  and  used  more  to  assist  the  boat  operator  and  divers 
in  station  keeping  during  servicing  than  as  a  required  locating  aid. 

10.  The  data  tapes  retrieved  from  the  gages  were  returned  to  WES  for 
processing.  The  wave  spectra  were  computed  via  a  discrete  Fast  Fourier  Trans¬ 
form  and  corrected  for  depth  attenuation  according  to  linear  wave  theory  as  a 
function  of  frequency.  The  pressure  response  factor  P  is  computed  as 

p  cosh  kh  . 

“  cosh  kD  u; 

where 

k  *  local  wave  number* 

h  ■  local  water  depth 

D  ■  height  of  the  wave  sensor  above  the  bottom 


i 


*  For  convenience,  symbols  are  defined  in  the  Notation  (Appendix  B) . 

t 

I 


The  spectral  energy  density  of  the  sea  surface  is  thus  related  to  the  pressure 
spectrum  by  the  equation 


V»  =  (gfr§)V>  <2> 

where 

subscripts  s,  p  =  surface  and  pressure  spectra,  respectively. 

The  significant  wave  height  H  is  computed  as 

mo 


1/2 


(3) 


where 

At  =  sampling  interval,  sec 
T  =  length  of  a  wave  record,  sec 

For  the  data  obtained  at  Saginaw  Bay,  T  =  1,024  sec;  for  the  data  obtained  in 
years  1979  and  1980,  At  =  0.5  sec;  for  1981,  At  =  1  sec.  The  sampling  inter 
val  was  increased  for  two  reasons.  First,  it  was  observed  in  the  wave  data 
obtained  during  1979  and  1980  that  there  was  little  energy  in  the  frequencies 
above  about  0.3  Hz;  second,  the  wave  climate  at  each  of  the  locations  could 
change  markedly  in  a  3-hr  period  which  was  the  length  of  time  between  wave 
records  acquired  during  1979  and  1980.  Thus  it  was  decided  tc  increase  the 
sample  interval  from  0.5  to  1.0  sec  and  decrease  the  interval  between  sample 
records  from  3  to  2  hr. 


PART  IV:  THEORY  OF  THE  SHALLOW-WATER  WAVE  MODEL 


11.  The  predictions  of  shallow-water  wave  characteristics  have  become  a 
focal  point  of  research  activities  across  the  world.  Because  construction, 
shipping,  and  dredging  operation  costs  have  drastically  increased  over  the 
years,  coastal  engineers  have  been  faced  with  more  accurately  defining  the 
shallow-water  wave  climate.  A  better  understanding  of  shallow-water  wave 
growth  and  transformation  mechanisms  is  slowly  evolving  through  controlled 
wave-measuring  programs  such  as  ARSLOE  (Vincent  and  Lichy  1981).  However,  not 
all  of  the  questions  have  been  answered,  and  it  will  take  some  time  before  all 
shallow-water  wave  transformation  mechanisms  are  quantified.  In  light  of  this, 
the  shallow-water  wave  modeling  technique  (SWWM)  employed  in  this  study  adopts 
"state-of-the-art"  mechanisms  currently  available.  The  main  intent  in  the 
development  of  the  SWWM  is  to  describe  the  physical  processes  as  accurately  as 
possible  while  simplifying  the  computational  procedures  to  a  degree  where 
shallow-water  wave  hindcasting  is  economically  feasible. 

12.  Hasselmann  et  al.  (1976)  introduced  a  parametric  model  of  wind-wave 
generation  relating  the  rate  of  energy  growth  to  nondimens ional  characteris¬ 
tics  of  the  wind  field.  The  energy  growth  (in  space  or  time)  is  governed  by  a 
self-similar  process  and  verified  through  extensive  prototype  data  (Hasselmann 
et  al.  1973,  1976).  In  these  studies,  the  dominant  energy  input  to  the  for¬ 
ward  face  of  the  spectrum  is  related  to  convergence  of  energy  flux  due  to  non¬ 
linear,  resonant  wave-wave  interactions  (Figure  7)  of  the  form  described  by 
Hasselmann  (1962).  Studies  by  Mitsuyasu  (1968,  1969)  and  Kitaigordskii  (1962) 
also  displayed  similar  results.  Although  these  formulations  were  developed 
for  deepwater  wave  conditions,  they  are  used  in  the  SWWM  because  the  only 
formulation  of  the  nonlinear  transfers  is  based  specifically  on  JONSWAP  type 
wave  spectra. 

13.  The  rate  of  wave  growth  under  ideal  conditions  of  fetch  limitations 
or  duration  limitation  and  a  stationary  wind  field  can  be  computed  (Hasselmann 
et  al.  1976).  For  growth  along  a  fetch  the  solution  is 

E  =  1.6  x  io"7  U2  -  (A) 

°  g 


and  for  growth  through  time,  it  becomes 
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Figure  7.  Schematic  representation  of  the  nonlinear 
wave-wave  interactions 


E  =  4.3  *  10-10  U,8/7  s'477  t10/7 


where  E  is  the  total  energy  resulting  from  a  wind  speed  U  (assumed  to  be 
o 

overwater  wind  conditions  adjusted  to  33-ft  elevation),  blowing  over  a  given 
fetch  length  F  .  The  gravitational  acceleration  is  denoted  by  g  ;  t  is 
the  time  since  the  wind  began  to  blow. 

14.  Two  additional  pieces  of  information  are  required  to  quantify 
the  distribution  of  Eq  given  in  the  form  of  an  energy  density  spectrum. 

The  nondimens ional  peak  frequency,  f  ,  and  the  Phillips'  equilibrium 


constant  a  (Phillips  1957)  are  shown  in  Figures  8  and  9.  These  parameters 

are  written  as 


a  =  0.0763T0,22 


(6) 


and 


2  =  3.5)f0,33 

m 


(7) 


where  $  is  the  nondimens ional  fetch  length 


$  =  Si 
2 
IT 


(8) 


15.  The  selection  of  fetch  (Equation  4)  or  duration  limited  conditions 
(Equation  5)  is  determined  from  the  following: 


t  .  =  5.37  x  io' 

min 


mf 


(9) 


where  is  the  minimum  duration  condition  and  Tg  is  the  significant 

wave  period  (U.  S.  Army  CERC  1981b)  given  by: 


T  =  7.54  -  tanh 
s  g 


^0.833 


tanh 


0.0379 


tanh 


(C 


0  833  ($t\ 


(10) 


where  h  is  the  mean  water  depth  along  F  . 

16.  If  t  ^  is  less  than  3  hr  (duration  of  each  input  wind  condition), 
then  Equation  5  will  be  used  to  compute  the  total  energy;  otherwise  Equation  4 
will  be  employed. 

17.  The  parameterization  of  the  wave  growth  is  restricted  such  that 
when  the  nondimens ional  peak  frequency  attains  a  value  of  0.13  or  less,  a 
fully  developed  sea  state  is  achieved  and  wave  growth  is  halted.  Over  long 
fetch  lengths  and  low  wind  speeds,  this  condition  can  occur  with  some  degree 
of  regularity.  Thus  Equations  4-8  are  then  redefined  by 


Figure  8.  Nondimensional  peak  frequency  f  as  a 
function  of  nondimensional  fetch  length  $ 
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Figure  9.  Phillips'  equilibrium  constant  a  ,  as  a 

function  of  nondimensional  fetch  length  $  ,  from 
Hasselman  et  al.  (1973) 


where  K  is  defined  as  the  nonvarying  parameters  (and  constants),  Q  is  de¬ 
fined  as  the  dependent  parameters,  and  is  recognized  as  the  independent 

parameters  (F  and  $)  found  in  Equations  4-8.  The  parameter  i  is  the  in¬ 
crement  counter.  After  each  discrete  fetch  length  F^  ,  the  nondimensional 
peak  frequency  is  evaluated  to  determine  if  ?  <  0.13  .  If  this  occurs  wave 

ID  “ 

growth  is  terminated,  and  wave  decay  is  initiated  for  the  remainder  of  the 
fetch  length.  Wave  decay  is  parameterized  following  the  work  conducted  by 
Bretschneider  (1952)  and  Mitsuyasu  and  Kimura  (1965)  for  f  the  peak  fre- 
quency  (where  f^  =  f^g/U)  while  the  total  energy  decay  rate  follows  that 
described  by  Jensen  (in  preparation) . 

18.  Wave  conditions  generated  in  Saginaw  Bay  also  must  consider  dis¬ 
persion  effects  resulting  from  finite  water  depth  conditions.  When  the  water 
depths  vary  from  F^  to  F^+^  »  the  conservative  transformation  mechanisms  of 
shoaling  and  refraction  must  be  considered.  Wave  shoaling  is  determined  from 
the  evaluation  of  group  speed  determined  by  linear  theory.  Wave  refraction  is 
neglected  under  the  assumption  that:  the  bottom  topography  is  assumed  to  be 
straight  and  parallel  for  every  fetch  length.  Considering  the  water  depths  in 
Saginaw  Bay  and  peak  wave  periods  (T  =  1/f^)  in  the  range  of  2  to  8  sec, 
wave-refraction  effects  (and  subsequent  "errors”)  would  be  on  the  order  of 

2  to  25  percent  as  shown  by  Figure  10,  designated  by  the  crosshatched  area. 
This  is  assuming  that  at  most,  the  angle  between  the  wave  crest  and  bottom 
contour  is  30  deg.  The  initial  direction  of  wave  propagation  is  limited  to 
18  angle  classes  at  10-deg  increments  (because  of  the  wind  data  employed  in 
this  study);  thus  the  accuracy  in  the  resultant  refracted  wave  condition,  by 
similarity,  also  would  be  constrained  to  the  18  angle  classes. 

19.  Finite  water  depth  conditions  also  lead  to  bottom  dissipation  ef¬ 
fects  on  the  growing  seas.  Energy  losses  associated  with  bottom  friction  are 
empirically  modeled  using  the  following  sets  of  equations  developed  by 
Bretschneider  and  Reid  (1954) : 

ff  E,  AF.  f4 
I  f  l  m 

K 

s 

where 

E  =  final  total  energy  at  F^ 

E.  =  original  total  energy  at  F._. 
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Figure  10.  Changes  in  wave  direction  due  to  refraction  on  slopes  with 
straight  and  parallel  depth  contours 


nondimens ional  friction  factor  (set  for  Saginaw  Bay  at  0.001) 
distance  of  wave  travel  within  the  discrete  fetch  length 


K  =  tanh 
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[2k.h.  T 
1  +  sinh  (2k.h.)J 


*  64n^  [  Ks  f 

"  3g2  L8inh  C2kihi)J 


k.  -  wave  number  (k.  =  =— ) 
x  i  L. 

x 

=  wavelength  evaluated  for  fQ 
h,  =  water  depth  at  F. 


20.  The  second  theoretical  aspect  of  SWWM  deals  primarily  with  the  dis¬ 
tribution  of  the  total  energy  (Eq)  in  the  form  of  a  one-dimensional  discrete 
frequency  spectrum  E(f^)  .  Through  the  use  of  similarity  principles,  Kitai- 
gordskii,  Krasitskii,  and  Zaslavaskii  (1975)  extended  Phillips’  deepwater  hy¬ 
pothesis  (Phillips  1958)  of  the  equilibrium  range  in  the  spectrum  of  wind¬ 
generated  surface  waves  to  finite  depth  conditions.  The  spectral  form  is 
defined  by 

E(f.)  =  ag2  (2n)'4  fT5  4>(u».  )  f.  >  f  (15) 

J  J  h  j  -  m 

where  E(f.)  is  the  energy  density  at  each  discrete  frequency  band,  f.  ,  and 
J  3 

♦(u^)  is  a  nondimens ional  function  dependent  on  given  by 


=  2nf j(h/g)' 


The  function  4»(uiji)  varies  from  1.0  in  deep  water  to  0.0  when  h  =  0.0  ,  as 
shown  by  Figure  11.  When  is  less  than  1.0,  ♦(u^)  can  be  approximated 

by: 


and  therefore, 


»<v=  K 


E(f.)  =  l  «gh  (2n)”2  f‘3 
J  ^  J 


f.  >  f 
j  -  m 


-5  -3 

or,  the  spectral  shape  changes  from  f  to  f  in  the  tail  of  the  energy 


density  spectrum,  and  more  importantly,  becomes  a  function  of  the  water  depth. 

21.  The  forward  face  of  the  spectrum  is  assumed  to  be  represented  by: 


E(f . )  =  org2  (2n)  ^  f  5  exp 
j  m 


-(f) 


* 


’(V 


f.  <  f 
J  ® 


(19) 


where  ♦' (ui^)  is  evaluated  from  the  defined  at  f^  .  Field  and  labo¬ 

ratory  data  by  Goda  (1974),  Thornton  (1977),  Ou  (1980),  Iwata  (1980),  and 
Vincent  (1981)  support  the  form  given  by  Equation  18.  The  verification  of 
Equation  19  can  be  found  in  Vincent*  and  is  supported  by  this  study. 

22.  The  parametric  representation  of  wave  growth  assumes  a  dynamic 
balance  between  atmospheric  sources  and  transfers  of  energy  resulting  from 
wave-wave  interactions  (Figure  6).  This  parameterization  was  based  on  deep¬ 
water  wave  conditions,  Hasselmann  et  al.  (1976).  During  this  study,  it  was 
determined  that  over  moderately  short  fetch  lengths  (10  to  20  n.m.),  this  deep 
water  growth  rate  expression  (Equations  4  and  5)  consistently  underpredicted 
the  total  energy  found  in  the  measured  data.  The  only  theoretically  consis¬ 
tent  location  to  add  the  energy  would  be  on  the  forward  face  of  the  spectrum 
(Figure  12).  The  function,  *s  t*ie  saturated  spectrum  based  on 

Equations  15  and  19,  and  E(f  ,h)ygjQ|jjgjj  is  the  spectrum  based  on  Eq  after 
wave  growth.  This  process  also  shifts  f^  to  a  lower  frequency  which  has 
been  noticed  in  field  data  (Vincent**) .  As  the  fetch  length  increases ,  the 
relative  amount  of  added  energy  decreases,  where  eventually  no  additional 
energy  is  incorporated  into  the  resulting  spectrum. 

23.  It  has  been  shown  that  the  water  depth  greatly  influences  the  spec 
tral  shape  and  in  so  doing  will  influence  the  maximum  wave  condition.  The 
parametric  formulation  follows  the  work  conducted  by  Vincent  (1981).  The 
depth  limiting  maximum  wave  condition  is  given  by, 


H  =  4 


(f)  df 


(20) 


*  Personal  communication,  C.  L.  Vincent  (1982a),  U.  S.  Army  CERC,  Fort 
Belvoir,  Va. 

**  Personal  communication,  C.  L.  Vincent  (1982b),  U.  S.  Army  CERC,  Fort 
Belvoir,  Va. 
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Figure  12.  Construction  of  the  final  energy  density  spectrum 
(solid  symbols)  caused  by  shallow-water  wave  generation 


where 


H  =  the  maximum  wave  condition 
m 

f  =  the  lower  frequency  bounding  the  total  energy  (equal  to  0.9  f^) 
E^Cf)  =  is  defined  from  Equation  12 

Integrating  Equation  20  the  absolute  limit  on  the  wave  condition  at  a  partic¬ 
ular  water  depth  is  obtained,  where 


c 


(21) 


24.  In  summary,  the  physical  process  governing  wave  generation  and 
transformations  has  been  theoretically  determined  using  available,  state-of- 
the-art  techniques.  It  must  be  emphasized  that  not  all  shallow-water  trans¬ 
formation  processes  have  (or  can  be)  measured  to  determine  their  relative 
effect  on  the  total  energy,  spectral  shape,  and  peak  frequency.  Therefore 
the  development  of  the  SWWM  as  employed  in  this  study  attempts  to  model  the 
physics  of  the  problem  in  a  general  sense  while  maximizing  computational 
efficiency. 


PART  V:  COMPARISONS 


25.  In  all  wave  hindcasting  studies,  comparisons  to  gage  measurements 
are  a  necessary  element  in  the  development  of  a  wave  model.  The  initial  cali¬ 
bration  test  was  conducted  on  a  data  set  that  contained  the  largest  wave  con¬ 
ditions  measured  during  the  wave  gaging  portion  of  this  study. 

26.  The  wave  gages  were  deployed  in  April  1981  shortly  after  the  bay 
became  ice-free.  Early  on  10  May  1981,  winds  began  to  increase  and  by  noon 
were  steady  at  25  to  30  mph  from  a  direction  of  about  40  deg  east  of  north. 

The  winds  held  a  remarkably  constant  speed  and  direction  for  about  a  day  and 

a  half  before  beginning  to  diminish.  The  predominant  direction  coincided  with 
the  axis  of  the  bay  and  alignment  of  the  gage  array,  the  most  favorable  condi¬ 
tion  for  generation  of  the  largest  waves  and  for  studying  changes  in  the  wave 
climate. 

27.  Because  of  the  constancy  in  the  wind  speed  and  direction  during  the 
period  10-12  May  1981,  any  variation  in  the  wave  climate  would  be  a  function 
only  of  the  wind  speed.  Figure  13  shows  the  wind  data  obtained  at  the  Saginaw 
Projects  Office  during  this  period  of  time.  The  anemometer  elevation  was 
about  60  ft  and  located  about  3  miles  south  of  the  disposal  facility. 

28.  Results  of  the  comparisons  for  gage  sites  1  and  2  are  shown  in 
Figures  14-17  where  H  is  the  significant  wave  height  and  T  is  the  peak 

m0  P 

period  defined  as: 


E(f)  df 


(22 


T 
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1.0 
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(23) 


There  is  a  slight  phase  difference  between  the  measured  and  hindcast  data 
sets.  This  is  due  to  the  SWWM  assumption  that  wave  conditions  are  generated 
instantaneously,  i.e.,  there  is  no  time-dependency  associated  with  the  effects 
of  wave  propagation  on  the  wave  climate.  The  small-scale  temporal  fluctua¬ 
tions  in  the  measured  data  cannot  be  simulated  in  the  SWWM  because  of  the  as¬ 
sumption  of  uniform  wind  conditions. 

29.  Energy  density  comparisons  are  also  made  between  the  measured  and 
hindcast  data.  It  is  helpful  to  recognize  that  a  spectral  representation  of 
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Figure  17.  Comparison  between  measured  and  hlndcast  peak 
wave  period  data  for  Station  2 

the  wave  phenomena  Is  an  estimate  of  the  actual  wave  conditions  existing  at  a 
specific  point  In  space  and  time.  To  Illustrate  this,  all  measured  wave  spec¬ 
tra  are  plotted  with  an  accompanying  90  percent  CHI-squared  distribution  con¬ 
fidence  band.  Although  numerous  comparisons  were  performed,  a  limited  number 
of  Station  1  results  are  presented.  The  five  energy  density  spectral  plots 
shown  in  Figures  18-22  represent  four  phases  in  the  10-12  Way  1981  storm;  Fig¬ 
ures  18  and  19  Illustrate  the  Initial  growth  phase.  Figure  20  illustrates  the 
initial  decay  phase.  Figure  21  illustrates  the  secondary  growth  phase,  and 
Figure  22  illustrates  the  final  decay  stage.  The  SWWM  spectral  results  are 
adjusted  in  time  to  compensate  for  the  lag  associated  with  the  wave  propaga¬ 
tion.  The  adjustment  varies  from  2  to  4  hr  depending  on  propagation  time. 

30.  The  energy  plots  (E(f)  versus  f)  are  plotted  in  a  nondimens lonal 
frequency  domain  defined  by  f/fm  where  fm  is  the  frequency  at  the  spectral 
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Figure  22.  Comparison  between  measured  and  hindcast  wave  spectra 
Station  1,  1000  hours,  11  May  1981 


peak.  The  greatest  discrepancy  in  the  comparison  of  measured  and  hindcast 
peak  frequencies  was  one  discrete  frequency  band  or  0.0156  Hz.  Figures  18-22 
show  that  from  initial  growth  to  final  decay,  the  SWWM  energy  density  spec¬ 
trum  follows  the  measured  spectrum  to  within  the  90  percent  confidence  limits. 
The  forward  face  of  the  measured  spectra  conforms  with  the  assumption  of  the 
spectral  shape  defined  in  Equation  19.  The  tail  of  the  hindcast  spectra  con¬ 
forms  with  the  measured  data  except  for  small-scale  temporal  variations  in  the 
measured  data  and  the  extreme  tail  of  the  distribution,  i.e.,  when  f/f  >  2  . 
The  discrepancy  between  the  measured  and  hindcast  spectra  when  f/f  >2  is 
a  characteristic  of  a  pressure  sensing  gage  system  such  as  used  in  this  study. 
Assuming  a  linear  transformation  exists  between  the  recorded  pressure  response 
and  the  free  surface,  individual  wave  component  frequencies  greater  than 
0.35  Hz  will  be  increasingly  damped  with  water  depth.  Thus  representation  of 
the  energy  density  at  frequencies  greater  than  =0.35  Hz  are  very  approximate. 
The  relative  energy  at  frequencies  greater  than  0.35  Hz  is  plotted  but  repre¬ 
sents  only  an  unknown  fraction  of  the  energy  in  a  particular  wave  record. 

31.  The  secondary  peaks  displayed  in  the  measured  spectra  may  be  caused 
by  several  different  mechanisms  such  as  secondary  wave  trains  propagating  into 
the  area  or  backscattering  due  to  variations  in  bottom  topography,  the  decompo 
sition  of  finite  amplitude  waves  into  secondary  harmonics,  or  an  artifact  of 
the  analysis  procedure. 

32.  Additional  comparisons  are  also  performed  on  three  storms.  For 
these  comparisons,  a  secondary  wind  information  data  source  was  needed  to 
supplement  the  Corps  Saginaw  Projects  Office  wind  data.  Wurtsmith  Air  Force 
Base,  located  near  Oscoda,  Michigan  (44.4  deg  N,  83.3  deg  W) ,  was  able  to  pro¬ 
vide  wind  data  for  the  additional  comparisons.  These  data,  along  with  the 
data  obtained  at  the  Corps  anemometer,  suggest  that  the  wind  conditions  over 
Saginaw  Bay  are  nearly  uniform.  Moreover,  the  Wurtsmith  wind  records  extend 
back  over  20  years,  thus  providing  a  long-term  data  base  suitable  for  hind¬ 
casting  extremes. 

33.  The  first  comparison  is  made  on  storm  conditions  during  4-7  May 
1981.  The  average  wind  speed  during  this  period  was  approximately  15  mph, 
with  maximum  winds  of  25  mph  (Figure  23).  The  wind  direction  is  shown  to 
vary  from  due  south  to  the  northeast.  Comparisons  between  the  Corps  and  the 
Oscoda  anemometers  demonstrate  that  the  wind  data  acquired  at  Wurtsmith  Air 
Force  Base  can  be  used  to  provide  an  accurate  representation  of  the  wind 


Figure  23.  Comparison  of  wind  records  obtained  at  Wurtsmith  AFB,  Oscoda, 
Michigan,  and  CE  Saginaw  Projects  Office,  Essexville,  Michigan 

climate  over  Saginaw  Bay.  Figure  22  shows  that  the  wind  data  acquired  at 
these  locations  are  virtually  the  same  after  adjustment  for  differences  in 

their  geographical  locations.  Figures  24  and  25  illustrate  the  results  of 

the  hindcast  and  measured  wave  conditions . 

34.  A  second  set  of  comparisons  is  made  for  a  period  of  moderate  winds 

which  occurred  during  8-10  October  1980.  The  wind  data  are  presented  in  Fig¬ 

ure  26,  and  again,  there  is  good  agreement  between  the  data  sets.  The  wind 
speed  during  this  period  averages  about  8  mph  with  a  maximum  of  about  15  mph. 
Results  of  the  hindcast  simulation  using  wind  data  from  the  Corps  anemometer 
are  presented  in  Figure  27.  The  relatively  large  discrepancies  are  probably 
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the  result  of  great  variability  in  the  wind  direction  during  periods  of  low 
wind  speeds . 

35.  Results  of  a  third  comparison  using  wind  data  shown  in  Figure  28 
are  shown  in  Figure  29.  The  wind  data  should  be  shifted  approximately  14  hr 
to  generate  a  wave  climate  to  coincide  with  that  observed.  The  necessity  for 
shifting  the  wind  data  is  probably  caused  by  an  incorrectly  annotated  wind 
record.  Once  adjusted  in  time  (not  shown  in  Figures  28  and  29),  the  discrepan 
cies  between  the  measured  and  hindcast  data  are  within  ±0.5  ft  for  the  signif¬ 
icant  height  and  +0.5  sec  for  the  peak  period. 

36.  Many  more  comparisons  were  performed  than  are  presented  here  with¬ 
out  adjustment  of  coefficients  or  modification  to  the  hindcast  model  to  ac¬ 
count  for  changes  in  location.  The  comparisons  convincingly  demonstrate  the 
ability  of  the  SWWM  to  accurately  hindcast  wave  conditions  in  Saginaw  Bay  to 
within  +0.5  ft  significant  height  and  +1.0  sec  peak  period.  Moreover,  the 
theoretically  derived  spectral  shapes  generally  conform  within  the  90  percent 
CHI-squared  confidence  limits  of  the  measured  spectra. 


DIRECTION.  DEG 


H/NDCAST  RESULTS 


PART  VI:  GENERATION  OF  DESIGN  WAVE  CONDITIONS 


37.  Specification  of  a  representative  wind  climate  is  necessary  for  the 
generation  of  design  wave  conditions  in  Saginaw  Bay.  The  largest  wind  waves 
in  the  bay  are  generated  by  synoptic-scale  systems  such  as  extratropical  cy¬ 
clones.  Smaller  scale  wind  systems  such  as  thunderstorms  and  squall  lines  are 
not  large  enough  and  do  not  persist  long  enough  to  generate  a  saturated  wave 
field.  Because  of  the  mesoscale  nature  of  the  wind  fields  of  interest,  one 

or  two  wind  records  on  the  periphery  of  the  bay  are  considered  adequate  to 
describe  the  wind  field  over  the  entire  bay. 

38.  A  recording  anemometer  system  was  installed  in  1980  at  the  CE  Sagi¬ 
naw  Projects  Office  located  near  the  mouth  of  the  Saginaw  Bay.  Informal  com¬ 
parisons  of  the  winds  measured  at  the  Saginaw  Projects  Office  with  those  mea¬ 
sured  on  the  bay  by  ships  operating  in  the  area  showed  good  agreement  after 
the  overwater  to  overland  transformation  was  made.  Records  obtained  at  the 
Saginaw  Projects  Office  were  also  compared  with  wind  records  obtained  at 
Wurtsmith  Air  Force  Base  located  near  Oscoda,  Michigan.  The  Wurtsmith  records 
are  particularly  useful  because  there  are  over  20  years  of  data  and  thus  pro¬ 
vide  a  basis  for  the  extremal  analysis.  Figure  23  shows  a  portion  of  a 
Wurtsmith  wind  record  during  a  period  in  which  data  from  the  CE  anemometer  at 
Essexville,  Michigan,  were  also  available.  The  direction  traces  show  good 
agreement  considering  the  distance  separating  the  stations.  The  speed  traces 
agree  less  well  but  part  of  the  disagreement  can  be  accounted  for  in  that  the 
CE  anemometer  at  Essexville  is  placed  at  a  considerably  higher  elevation  than 
the  Wurtsmith  anemometer. 

39.  Thus  while  the  CE  anemometer  at  Essexville  provided  the  wind  data 
for  verification  of  the  SWWM,  the  wind  records  obtained  at  Wurtsmith  provide 
the  basis  for  the  extremal  analysis. 

40.  In  PART  V,  the  SWWM  was  verified  and  a  capability  to  hindcast  the 
time-history  of  wave  activity  during  individual  storms  was  demonstrated.  This 
capability  permits  the  computation  of  extremal  statistics  if  sufficient  wind 
data  are  available.  One  of  the  major  tasks  of  this  study  was  locating  suf¬ 
ficient  wind  data  to  confidently  perform  the  analysis. 

41.  In  PART  V  it  was  demonstrated  that  wind  data  obtained  at  Wurtsmith 
could  be  used  to  reliably  hindcast  wave  conditions  in  Saginaw  Bay.  Twenty 
years  of  wind  data  (1950-1970)  were  requested  and  obtained  from  Wurtsmith. 


The  extremal  hlndcasts  were  performed  under  the  following  assumptions: 

Wind  speed  and  direction  are  held  constant  for  each  event. 

1>.  The  waves  are  generated  under  fetch-limited  conditions. 

The  first  step  in  the  extremal  analysis  was  to  determine  the  wave  heights  pro¬ 
duced  by  constant  wind  speed  conditions  for  each  azimuthal  angle  band.  The 
angle  bands  selected  cover  10  deg  of  azimuth.  Figure  30  shows  the  results  of 
this  analysis;  the  highest  Hm  values  are  generated  by  winds  from  the  20-  to 
40-deg  azimuthal  angle  bands.  ^The  absolute  maximum  wave  condition  (Vincent 
1981,  Equation  21)  of  7.9-ft  significant  wave  height  is  generated  when  the 
wind  speed  attains  50  mph  in  these  direction  bands.  There  is  a  significant 
drop  in  the  wave  height  in  the  320-deg  angle  band  caused  by  the  limited  fetch 
length  and  water  depths  that  exist  for  that  angle  class.  The  data  presented 
in  Figure  30  should  not  be  used  for  design  purposes.  The  Hm  results  are 
generated  for  very  specific  wind  conditions  that  may  not  have^occurred  in 
nature.  The  curves  are  used  in  this  case  to  further  quantify  extreme  storm 
conditions  for  given  wind  speeds  and  principal  wind  directions.  Once  the 
maximum  wave-producing  angle  bands  are  known,  it  becomes  a  matter  of  search¬ 
ing  the  20-year  wind  record  for  wind  conditions  in  the  20-  to  40-deg  angle 
bands  and  wind  speeds  greater  than  25  knots.*  Of  a  total  of  153,298  observa¬ 
tions,  631  recorded  wind  speeds  greater  than  25  knots  for  all  direction  bands 
(Table  1) .  Of  the  631  observations,  five  wind  conditions  were  selected  that 
as  a  group,  produced  the  most  severe  wave  conditions.  The  "most  severe"  is 
not  to  say  that  no  greater  conditions  existed;  simply  that  the  anemometer  may 
not  have  been  in  continuous  operation  during  such  periods.  The  wind  condi¬ 
tions  producing  the  most  severe  wave  conditions  are  shown  in  Table  2. 

42.  The  time-histories  of  the  storms  are  provided  in  Appendix  A  (Fig¬ 
ures  A1-A10).  The  wind  speed  and  direction  were  averaged  over  2-hr  intervals. 
Wind  data  for  the  five  storm  conditions  were  used  to  drive  the  SWWM  to  compute 

H  and  T  .  The  winds  were  corrected  for  differences  in  anemometer  heights 
mQ  P 

ana  converted  to  overwater  winds  according  to  the  procedure  described  by 
Resio  and  Vincent  (1976)  and  U.  S.  Army  CERC,  CETN-I-5  (1981a).  The  winds 
and  corresponding  wave  conditions  are  displayed  in  Figures  A1-A10.  Table  3 
is  a  tabulation  of  the  maximum  significant  wave-height  conditions  computed 
using  the  SWWM. 


*  Wurtsmith  Air  Force  Base  recorded  the  wind  data  in  units  of  whole  knots. 
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Table  1 

Percent  Occurrence  of  Wind  Conditions  for  Wurtsmlth  Air  Force  Base  1950-1970 

(Percent  Occurrence  Times  1000) 


Note:  Total  number  of  observations  greater  or  equal  to  25  knots  ■  631. 
Total  number  of  observations  from  1950-1970  ■  153,298. 


Table  2 

Storm  Conditions  for  Extreme  Analysis 


43.  Although  arguments  can  be  made  for  ranking  the  storms  by  either 
or  H  ,  the  latter  ranking  is  used  in  this  study  because  of  the  exagger¬ 
ated  influence  of  possible  wind  errors.  Consider  the  maximum  H  ,  7.1  ft, 

m0 

found  in  the  first  storm  vrtiich  was  generated  using  a  windspeed  of  37  knots. 

The  37-knot  wind  speed  was  calculated  by  averaging  wind  speeds  +1.0  hr  around 
a  given  time  t  as 


ut  ‘  t(Vi  ♦  °t  *  Vi) 


The  three  values  employed  for  ut_^  »  ut  »  and  Ut+1  were  25,  50,  and  36  knots 
respectively.  It  appears  that  the  50-knot  wind  speed,  of  which  there  is  only 
one  value  in  the  20-year  record,  may  be  incorrect  or  a  short  gust  unrepresen¬ 
tative  of  the  synoptic  scale  winds.  Thus  Table  4  shows  the  storms  ranked 
according  to  the  H  value  with  return  period 

mo 


1  +  NT 


where 


Tr  *  the  return  period 
R  -  the  storm  rank 

NY  *»  the  duration  of  the  data  base  in  years 
These  data,  along  with  the  return  period,  are  displayed  in  Figure  31. 


Table  4 

Storm  Ranking  and  Return  Period 


Storm  Number 


Return  Period,  year 


44.  TVo  additional  storm-related  values  are  presented  in  Fig¬ 

ure  31:  the  maximum  wave  condition  occurring  during  the  10  May  1981  storm 


5 

SIGNIFICANT  HEIGHT,  FT 


6 


(measured),  and  the  design  maximum  wave  condition  (NCE*)  used  for  the  construe 

tion  of  the  retaining  structure.  It  is  evident  that  both  H'  conditions 

m0 

are  considerably  smaller  than  the  H  conditions  hindcast  for  maximum  storm 

m0 

related  H  data.  Comparison  of  the  two  data  points  versus  the  mean  H 

m0  m0 
data  shows  that  the  resulting  return  period  would  be  approximately  6  years 

(10  May  1981  storm)  and  15  years  (design  wave  condition).  Although  the  10  May 

1981  storm  generated  the  largest  H  conditions  recorded,  they  are,  in  com- 

m0 

parison  with  historical  data,  of  relatively  small  magnitude. 

45.  Additional  estimates  of  extreme  wave  conditions  can  be  obtained 
employing  an  assumed  statistical  distribution  of  the  waves  within  a  given 


recording  period.  A  Rayleigh  distribution  has  normally  been  assumed  to  rep¬ 
resent  the  observed  wave  conditions  (U.  S.  Army  CERC  1977).  From  this  as¬ 


sumption,  a  relationship  is  obtained  relating  Hm  to 
is  the  average  of  the  highest  1  percent  of  all  waves),  or 


(where  Hj/1()& 


H,,,..  =  eH 
1/100  mQ 


where  e  is  equal  to  1.67.  Care  must  be  employed  when  using  the  above 

equation  since  the  H  results  are  theoretically  limited  (and  physically 

®0 

limited)  by  the  waterdepth  (Equation  21)  and  also  saturation  conditions 
(Figure  31).  Therefore  will  be  limited  in  a  similar  manner  as 

When  H  approaches  its  theoretical  maximum  condition  (approximately  u 
7.9  ft)  9he  assumed  Rayleigh  distribution  will  no  longer  hold;  the  tail  of 
the  distribution  would  be  truncated  at  some  wave  height  (Baitjes  1972  and 
Goda  1975).  When  this  occurs  e  will  diverge  from  the  assumed  wave  of  1.67 
and  approach  1.0. 


*  P.  McCallister.  1979.  Letter,  "Subject:  Saginaw  Diked  Disposal  Facility, 
U.  S.  Army  NCE,  Detroit,  Mich. 


PART  VII:  SUMMARY  AND  CONCLUSIONS 


46.  A  wave-measurement  system  was  deployed  in  Saginaw  Bay,  Michigan,  in 
the  fall  of  1979  and  operated  during  the  ice-free  parts  of  the  spring  and  fall 
seasons  of  1980  and  1981.  Concurrently,  development  of  a  numerical  wave  hind- 
cast  model  capable  of  simulating  wind-wave  generation  in  shallow  water  was 
completed . 

47.  The  most  severe  and  definitive  storm  during  the  measurement  period 
occurred  in  May  1981.  While  the  wind  speeds  recorded  during  this  storm  were 
not  particularly  high  (30  to  35  mph) ,  they  were  very  uniform  in  terms  of  both 
speed  and  direction.  Moreover,  the  predominant  direction  was  parallel  to  the 
major  axis  of  the  bay  thereby  allowing  the  greatest  fetch  length.  The  condi¬ 
tions  that  occurred  during  this  storm  resulted  in  a  near  ideal  wave  data  set 
for  verifying  the  Shallow-Water  Wave  Model  (SWWM) .  All  three  gages  performed 
almost  perfectly  and  data  recovery  was  nearly  100  percent.  The  near  station¬ 
ary  wind  field  permitted  a  definitive  test  of  correct  simulation  of  growth 
with  fetch  and  growth  with  time  of  the  SWWM.  Results  of  the  verification  are 
presented  herein  and  demonstrate  the  ability  of  the  SWWM  to  reproduce  both  the 
total  energy  and  peak  frequency  of  the  measured  spectra. 

48.  Review  of  20  years  of  wind  data  obtained  at  Wurtsmith  Air  Force 

Base  near  Oscoda,  Michigan,  disclosed  five  storms  that  occurred  between  1950 

and  1970  suitable  to  use  as  a  basis  for  an  extremal  analysis.  In  the  20  years 

of  recorded  data,  there  was  one  occurrence  of  a  50-mph  wind  speed.  Hindcast 

simulation  using  this  wind  speed  showed  the  wave  field  to  reach  a  saturated 
condition  at  about  7.9-ft  maximum  wave  height.  Thus  it  appears  the  maximum 
wave  height  to  be  expected  at  the  disposal  facility  is  about  8  ft. 

49.  There  have  been  informal  reports  of  sustained  wind  speeds  greater 
than  60  mph  measured  prior  to  the  inception  of  the  wave-measurement  program. 

It  is  recommended  that  the  wind  measurement  system  installed  at  the  Saginaw 
Projects  Office  continue  in  operation  indefinitely.  Wind  data  of  this  nature 
will  provide  very  valuable  information  on  the  frequency  of  occurrence  of  very 
high  wind  speeds  and  may  lead  to  revision  of  the  table  of  return  periods  for 
large  wave  heights. 
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Figure  A2.  Hlndcast  significant  wave  height  and  peak  wave  period 
for  atom  of  30  October  to  1  November  1951 
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Figure  A3.  Wind  data  used  for  hlndcast  of  storm  of 
18-20  November  1958 
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Hindcast  significant  wave  height  and  peak  wave  period 
for  storm  of  25-27  November  1958 
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Figure  A10.  Hindcast  significant  wave  height  and  peak  wave  period 

for  storm  of  8-10  April  1968 
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Wave  number  k.  =  — 
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Nonvarying  parameters 


Wavelength  evaluated  for  f 


m 


Duration  if  data  base  in  years 
Pressure  response  factor 
Dependent  parameters 
Storm  rank 

Time  since  wind  began  to  blow 
Minimum  duration  condition 


T  Length  of  wave  record 
Peak  wave  period 
Return  period 

Ts  Significant  wave  period 
U  Wind  speed 

$  Nondiaensional  fetch  length 
a  Phillips 1  equilibrium  constant 
Afc  Sampling  interval 

AF\  Distance  of  wave  travel  within  discrete  fetch  length 
Independent  parameters 

2 

a  Variance  about  the  mean  H  (or  H  ) 

"°  4  \ 

4-  Friction  function,  units  of  T  -  L 
4(ui^)  Nondimens ional  function  dependent  on  u)^ 


Subscripts 


Surface  and  pressure  spectra,  respectively 


